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Abstract 

We consider a class of Hamiltonians with three degrees of freedom that can be 
mapped into quasi-periodically driven pendulums. The purpose of this paper is to 
determine the threshold of the break-up of invariant tori with a specific frequency 
vector. We apply two techniques : the frequency map analysis and renormalization- 
group methods. The renormalization transformation acting on a Hamiltonian is a 
canonical change of coordinates which is a combination of a partial elimination of 
the irrelevant modes of the Hamiltonian and a rescaling of phase space around the 
considered torus. We give numerical evidence that the critical coupling at which 
the renormalization transformation starts to diverge is the same as the value given 
by the frequency map analysis for the break-up of invariant tori. Furthermore, we 
obtain by these methods numerical values of the threshold of the break-up of the 
last invariant torus. 

Key words: Invariant tori, Renormalization, Hamiltonian systems 
PACS: 05.45.Ac, 05.10.Cc, 45.20.Jj 



Preprint submitted to Physica D 



8 February 2008 



1 Introduction 



For Hamiltonian systems, the persistence of invariant tori influences the global 
properties of the dynamics. The study of the break-up of invariant tori is 
thus an important issue to understand the onset of chaos. For two degrees 
of freedom, there are several numerical methods to determine the threshold 
of the break-up of invariant tori : for instance, Greene's criterion [1], ob- 
struction method [2], converse KAM [3,4], frequency map analysis [5-8], or 
renormalization-group methods [9-11]. 

In this article, we propose to compute this threshold for a one-parameter fam- 
ily of Hamiltonians with three degrees of freedom and for a specific frequency 
vector, by two techniques : by frequency map analysis and by renormaliza- 
tion. The frequency map analysis is valid for any dimension, and has been 
applied to systems with a large number of degrees of freedom [5]. The set-up 
of renormalization-group transformations is also possible for any dimensions 
in the framework of Ref. [12], but only systems with two degrees of freedom 
have been investigated numerically. 

We describe the renormalization-group transformation and we implement it 
numerically for the spiral mean torus. The result is that the values of the crit- 
ical coupling given by the renormalization coincide up to numerical precision 
with the thresholds of the break-up of the spiral mean torus (of dimension 
3) given by frequency map analysis. The two methods we compare are com- 
pletely independent, both conceptually and in their practical realizations. The 
frequency map analysis is based on the analysis of trajectories, while the renor- 
malization is based on a criterion of convergence of a sequence of canonical 
t r ansf or mat ions . 

We conjecture, on the basis of this numerical result, that the renormalization- 
group transformation converges up to the critical surface (the set of Hamilto- 
nians where the torus of the given frequency is critical, i.e. at the threshold 
of its break-up), at least in a region of the critical surface of the Hamiltonian 
space where critical couplings are small enough (in order that the elimination 
procedure is well-defined [13]). 

We consider a class of Hamiltonians with three degrees of freedom written in 
terms of actions A = (Ax, A 2 , A 3 ) E M 3 and angles cp = (</?i, (p 2 , ^3) £ T 3 (the 
3-dimensional torus parametrized, e.g., by [0, 27r] 3 ) 

H e (A,<p)=H (A)+eV(<p), (1) 

where e denotes the coupling parameter. In this article, we consider the par- 
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ticular class of models for which the integrable part H is given by 

H (A) = Lj -A + ^(n-A) 2 , (2) 

where oj is the frequency vector of the considered invariant torus, and Q is 
another constant vector non-parallel to u> . We suppose that u> is incommen- 
surate, i.e. there is no nonzero integer vector v such that u> ■ v = 0. 
Since the quantity ll>q ■ cp is conserved (where uJq denotes a vector orthogonal 
to Q and to u; ), one can show (even if uJq ■ is not a function on the three- 
dimensional torus) that this model (l)-(2) is intermediate between two and 
three degrees of freedom; in appropriate coordinates it can be interpreted as 
one degree of freedom driven by a multi-periodic force with incommensurate 
frequencies u> . In particular, invariant tori in this intermediate model act as 
barriers in phase space (limiting the diffusion of trajectories) in a similar way 
as for two degrees of freedom Hamiltonian systems. We analyze in this article 
the break-up of invariant tori with spiral mean frequencies for this particular 
type of models, by choosing a special form of the perturbation (see section 3), 
such that the model is equivalent to a pendulum driven by two periodic forces 
with incommensurate frequencies. The method is however applicable to any 
perturbation and to the case of full three degrees of freedom [12,19]. 
We are interested in the stability of the torus with frequency vector u> . For 
the unperturbed Hamiltonian H Q , this torus is located at A = 0. Kolmogorov- 
Arnold-Moser (KAM) theorems were proven for Hamiltonians (1) provided 
that ojq satisfies a Diophantine condition [14]. This theorem shows the exis- 
tence of the torus with frequency vector ujq for a sufficiently small and smooth 
perturbation eV. The invariant torus is a small deformation of the unperturbed 
one. The existence of the torus outside the perturbative regime is still an open 
question even if efforts have been made to increase lower bounds for specific 
models (for a two dimensional model, see Ref. [15,16]). Conversely, for suf- 
ficiently large values of the coupling parameter, it has been shown that the 
torus does no longer exist [3,4]. The aim of this paper is to determine e c such 
that H £ has a smooth invariant torus of the given frequency for |e| < e c , and 
does not have this invariant torus for \e\ > e c . 

The invariant torus we study (named the spiral mean torus) has the frequency 
vector 

u> = (a 2 , a, 1), 

where a is the spiral mean, i.e. the real root of a 3 = a + 1 (am 1.3247). From 
some of its properties, a plays a similar role as the golden mean in the two 
degrees of freedom case [17]. The analogy comes from the fact that one can 
generate rational approximants by iterating a single unimodular matrix N. In 
what follows, we call resonance an element of the sequence {v k = N k ~ 1 Ui, k > 
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1} where V\ = (1,0,0) and 



N 



The word resonance refers to the fact that the small denominators u> ■ u k 
appearing in the perturbation series or in the KAM iteration, tend to zero 
geometrically as k increases (u>o • v k = a 3 ~ k — > as k — > oo). We notice that 
Wo is an eigenvector of N, where N denotes the transposed matrix of N. One 
can prove [12] that u> satisfies a Diophantine condition of the form : 

\(jl> ■ v\ > C\l>\~ 2 , 

where \u\ = (\vi\ 2 + \v 2 \ 2 + l^j 2 ) 1 / 2 , and c « 0.6. 



2 Renormalization-group transformation 

The renormalization transformations are defined for a fixed frequency vector 
ujq, and contain a partial elimination of the irrelevant modes (the non-resonant 
part) of the perturbation, and a rescaling of phase space. The elimination of 
irrelevant modes is performed by iterating a change of coordinates as in KAM 
theory. We remark that other perturbative techniques can be used instead, 
leading to similar results. The rescaling of phase space combines a shift of 
the resonances, a rescaling of time, and a rescaling of the actions. The aim is 
to change the scale of the actions to a smaller one (and a longer time scale). 
This renormalization can be thought as a microscope in phase space. The non- 
resonant modes are the ones which affect the motion at short time scales, and 
can be dealt with averaging methods. We define the non-resonant modes to 
be the ones that satisfy the inequality 

1 , , 

This set of modes, denoted I~, is the interior of a cone around the Wo-direction 
in the space of 3-dimensional vectors, with angle tc/A. We define the resonant 
modes as the Fourier modes which do not satisfy the condition (3), i.e. this set, 
denoted I + , is the complement of I in Z 3 . Since v k does not satisfy Eq. (3) 
for k > 1, I + contains the resonances that produce small denominators in the 
perturbation series or in the KAM theory. The "frequency cut-off" (between 
resonant and non-resonant modes) restricts the Fourier modes that can be 
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eliminated in one renormalization step, without running into small denomina- 
tor problems (the non-resonant modes). As it is common with cut-offs, there is 
not a single "natural" choice. More generally, other choices in the splitting of 
{e tl/ <p } into resonant and non-resonant modes should lead to the same results 
provided, e.g., that the ratio |i/|/|u;o • v\ is bounded on J", and that the shift 
of the resonances contracts non-zero vectors v in I + and maps them into I 
after a finite number of iterations of the transformation [11,12]. 



The transformation, acting on a Hamiltonian H of the form 

H(A,<p) = H (A) + V(n-A,<p), (4) 
where H is given by Eq. (2), combines four steps: 

(1) We shift the resonances v k+1 h- > v k : We require that the new angles <p' 
satisfy 

cos(V fe+ i • (f) = cos(v k ■ </?'), 
for k > 1. This is performed by the linear canonical transformation 

(A,cp)^(A',cp') = (N- 1 A,Ncp). 

Since N is an integer matrix with determinant one, this transformation pre- 
serves the T 3 -structure of the angles. We notice that the resonance v\ is 
changed into u = (0,1,1) which satisfies the condition (3), i.e. it is a non- 
resonant mode : Some of the resonant modes are turned into non-resonant 
ones by this linear transformation. 

This step changes the frequency a» into Nu> = o" _1 ^o (since uj is an eigen- 
vector of N by construction), and the vector f2 into Nfl. In order to keep a 
unit norm, we define the image of Q by 

fl' = . 5 

||iVO|| V 1 

(2) We rescale the energy (or equivalently time) by a factor a (i.e. we multiply 
the Hamiltonian by a), in order to keep the frequency fixed at ujq. 

(3) We rescale the actions : 

H'(A,tp) = \H ^j, 

such that the mean- value of the coefficient of the quadratic term in H' is equal 
to (Q' ■ A) 2 /2. This normalization condition is essential for the convergence of 
the transformation. After Steps 1, 2 and 3, the Hamiltonian expressed in the 
new variables is 

H'( A, <p) = XaH (jNA, N~ V) ■ (6) 
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For H given by Eq. (4), this expression becomes 

H\A, (p) = oj ■ A + ^||jVfi|| 2 (fi' ■ A) 2 + XaV (M^fl' ■ A, N' 1 ^ . 

(7) 

Thus the choice of the rescaling in the actions (Step 3) is 

\ = (r\\Nfl\\ 2 (l + 2(V^)), (8) 

where (V^) denotes the coefficient of the mean- value of the quadratic part, 
in the representation : 

V(Qf.A,ip) = ^V^(tp)(sr.A)'. 

i>o 

(4) We perform a near-identity canonical transformation Uh that eliminates 
completely the non-resonant part of the perturbation in H' . This transforma- 
tion satisfies the following equation : 

r(H'oU H ) = 0, (9) 

where I~ denotes the projection operator on the non-resonant modes acting 
on a Hamiltonian H as 

I-H(A,cp)= £ H v {Ay u(p - (10) 

Equation (9) is solved by a Newton method. We iterate a change of coordinates 
as in KAM theory that reduces the non-resonant modes of the perturbation 
from e to e 2 . One step of the elimination is performed by a Lie transformation 
Us : (A, if) i— > (A', if'), generated by a function S(A, cp). The expression of a 
Hamiltonian H in the new coordinates is given by 

HoU s = exp(S)H = H + {S,H} + ±{S, {S, #}} + •••, (11) 

where { , } is the Poisson bracket between two scalar functions of the actions 
and angles: 

it -i df dg df dg 

{f ' 9} = ^p-dA-dA-^ (12) 

and the operator S is defined as SH = {S,H}. The generating function S 
is chosen such that the order e of the non-resonant part of the perturbation 
vanishes. We construct recursively a series of Hamiltonians H n , starting with 
Hi = H', such that the limit Hoo is canonically conjugate to H' but does 
not contain non-resonant modes, i.e. I'H^ = 0. One step of this elimination 
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procedure, H n \— > H n+1 , is done by applying a change of coordinates U n such 
that the order of the non- resonant modes of H n+ i = H n o U n is e^, where 
e n denotes the order of the non-resonant modes of H n . At the n-th step, the 
order of the non-resonant modes of H n is e^" , where e is the order of the 
non-resonant modes of H' . If this procedure converges, it defines a canonical 
transformation Uh — U\ o U2 o ■ ■ ■ o U n o • • • , such that the final Hamiltonian 
Hoc — H' o Uh does not contain any non-resonant mode. 
The specific implementation of one step of this elimination procedure will be 
described more explicitly in the next sections. We will discuss two versions of 
this transformation : The first one is a renormalization for Hamiltonians in 
power series in the actions, and the second one is a slightly different version 
which eliminates only the non-resonant modes of the constant and linear part 
in the actions of the rescaled Hamiltonian H', and which allows us to define 
a renormalization within a space of quadratic Hamiltonians in the actions, 
following Thirring's approach of the KAM theorem [18]. 
It has been proven in Ref. [12] that, for a sufficiently small non-resonant part 
of the perturbation, the transformation Uh is a well-defined canonical trans- 
formation such that a Hamiltonian expressed in these new coordinates does 
not have non-resonant modes. The domain of definition of Uh has been ex- 
tended to some non-perturbative domain in Ref. [13]. 

Concerning the quadratic case, we lack at this moment a theoretical back- 
ground to prove an analogous theorem. The convergence of the elimination 
procedure in both cases outside the perturbative regime (e small) is observed 
numerically. 

In summary, the renormalizat ion-group transformations we define act as fol- 
lows : First, some of the resonant modes of the perturbation are turned into 
non-resonant modes by a frequency shift and a rescaling. Then, a KAM-type 
iteration eliminates these non-resonant modes, while slightly changing the res- 
onant modes. 



2. 1 Renormalization scheme for Hamiltonians in power series in the actions 

We define in this section, one step of the elimination of the non-resonant modes 
of the rescaled Hamiltonian H' . The renormalization transformation acts on 
the following family of Hamiltonians 

00 

H(A, <p) = uj ■ A + £ V®(<p)(n ■ AY, (13) 

j=0 

with = 0. We suppose that (V^) is nonzero (in order to have a twist 

direction in the actions). The approximations involved in this transformation 
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are of two types : we truncate the Fourier series of the functions i.e. we 
approximate a scalar function / of the angles by 



f [ - L] (v) = E (14) 

vec L 

where Cl = \y = (v\, ^2, ^3) G ^ 3 |maxj \vi\ < L}, and we also neglect all the 
terms of order O ((fi • A) J+1 ^j that are produced by the transformation. 

One step of the elimination procedure H 1— > H = H o is constructed as 
follows : We consider that depends on a small parameter e, such that 
is of order e. We define H , the integrable part of H as 

tfo(A)=u> -A+(V( 2 )>(ft-A) 2 , (15) 

and the perturbation of H is denoted V' = H-H = V- (V {2) )(Q ■ A) 2 . In 
order to eliminate the non-resonant modes of to the first order in e, we 
perform a Lie transformation U s : (A, 1— > (A', y>') generated by a function 
5 of order e and of the form 

j 

v?) = i E V 0) (^)(O • A) J ' + aQ, ■ cp. (16) 

The first terms of H are H + V + {5, # } + {5, V'} + 0(e 2 ). The function 5 
is determined by imposing that the order e vanishes : 

r{s,H } + tv + r{s,i + v'} = o. (17) 

This condition determines the non-resonant modes of S. For the resonant ones, 
we choose I + S = 0. Thus we notice that the mean value of {S, I + V'} is zero. 
The constant a eliminates the linear term in the (ft ■ A)-variable, (V^), by 
requiring that ({S, H }) + (V^)fl ■ A = 0: 

(V<V) 

(18) 



2{1 2 (VW)' 

We solve Eq. (17) by a Newton method with an initial condition which satisfies 
I~{iS f , H } + I~V = 0, since I + V is expected in general to be small. Then we 
compute H = H o Us by calculating recursively the Poisson brackets S k H = 
SS^H, for k > 1. Denoting H k = S k H, H becomes 

OO TT 

H = T.% (19) 



8 



2.2 Thirring's scheme for quadratic Hamiltonians 

We consider the following family of quadratic Hamiltonians in the actions and 
described by three scalar functions of the angles 

H(A, if) = u)q ■ A + m{<p)(n ■ Af + g(<p)Q ■ A + f(<p), (20) 

The KAM transformations are constructed such that the iteration stays within 
the space of Hamiltonians quadratic in the actions [18]. In order to prove 
the existence of a torus with frequency vector oj for Hamiltonian systems 
described by Eq. (20), it is not necessary to eliminate m, but only g and 
/ (the main point is that the torus with frequency vector u>o is located at 
A = for any H with / = g = 0, even if H is not globally integrable). The 
elimination of / and g can be achieved with canonical transformations with 
generating functions that are linear in the action variables, and thus map the 
family of Hamiltonians (20) into itself. This is very convenient numerically, as 
one only works with three scalar functions m, g and /. The only approximation 
involved in the numerical implementation of the transformation is a truncation 
of the Fourier series of these functions, according to Eq. (14). 
In this section, we describe one step of the elimination of the non-resonant 
modes H \— > H = H o Us- We assume that g and / depend on a (small) 
parameter e, in such a way that I~g and I~f are of order e. The idea is to 
eliminate the non-resonant modes of g and / to first order in e, at the expense 
of adding terms that are of order 0(e) in the resonant modes and of order 
0(e 2 ) in the non-resonant modes. This is performed by a Lie transformation 
Us : (A, (p) i— > (A', ip') generated by a function S of order e linear in the 
action variables, of the form 

S(A,<p) = Y((p)n-A + Z((p)+aQ-(fi , (21) 

characterized by two scalar functions Y, Z, and a constant a. The expression 
of the Hamiltonian in the new variables is obtained by Eq. (11). A consequence 
of the linearity of S in A is that the Hamiltonian H is again quadratic in the 
actions, and of the form 

H(A, <p) = w • A + m(cp)(n ■ A) 2 + g(<p)Q ■ A + f(cp). (22) 

This can be seen by this simple argument : Given a quadratic function in the 
(f2 • A)-variable and a linear function S, then SH is again quadratic in the 
(ft ■ A)-variable. The derivatives dH/dA and dS/dtp are linear in the (ft ■ A)- 
variable, and dH/d<p is quadratic, while dS/dA is constant in this variable. 
Therefore, SH given by Eq. (12) is quadratic. Consequently, iterating this 
argument, exp(S)H is also quadratic. We notice that the vector f2 remains 
unchanged during each step of the elimination. 

The functions Y, Z, and the constant a are chosen in such a way that 
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and I~f vanish to order e. The constant a corresponds to a translation in the 
actions, which has the purpose of eliminating the mean value of the linear 
term in the (ft ■ A)-variable. Then we express H by calculating recursively the 
Poisson brackets S k H = SS 1 * -1 !!, for k > 1, like in the previous section. 
We notice that for some purposes it is more convenient to eliminate also the 
non-resonant part of m, together with the one of g and / (see the remarks in 
Refs [12,11]). But such elimination procedure generates arbitrary orders in the 
Q ■ A-variable and this leads to the first version of the transformation. The 
advantage to work with this second version of the elimination procedure is 
that the Hamiltonians (20) are described by only three scalar functions of the 
angles. Thus it is numerically more efficient since the renormalization map is 
of lower dimension. Concerning the renormalization transformation for Hamil- 
tonians in power series in the actions, we truncate the Fourier series of each 
scalar function of the angles with a cut-off parameter L, and the Taylor series 
with a cut-off parameter J. For fixed L and J, the dimension of the renormal- 
ization map is equal to (J + 1)(2L + l) 3 + 2. Concerning the renormalization 
defined for quadratic Hamiltonians, the renormalization map is of dimension 
3(2L + l) 3 + 2. 

Remark : Another advantage to work the Thirring's version of the renormaliza- 
tion is that it can be generalized more easily to non-degenerate Hamiltonians 
H(A, if) = H (A) + V(A, ip) with Rank d 2 H /dA 2 > 1. (see Ref. [19]). 



3 Determination of the critical coupling 

We consider the following quasi-periodically driven pendulum model : 

H = \r? — p + e (cos x + cos (re + fit) + fi cos(x + z/ 2 t)) , (23) 

with the frequencies V\ = a + 1, u 2 = o 2 + 1, and /i a real parameter. This 
model can be mapped into the following degenerate Hamiltonian system with 
three degrees of freedom : 

H £ (A, ip) = u ■ A + • A) 2 + ef(<p) , (24) 
where fl = (1,1,-1) and the perturbation / is given by 

f((fi) = At cos y?i + cos y? 2 + cos tp s , (25) 
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This can be seen by considering the three angles ip[ = x, tp' 2 = v\t mod 2-7T, 
and ip 3 = v 2 t mod 2n. The Hamiltonian (23) becomes : 



H = -A' 2 - A[ + v x A 2 + v 2 A' 3 + e (cos (p[ + cos((p[ + (p 2 ) + n cos(v? / 1 + <// 3 )) , 

(26) 

where we added v\A' 2 + v 2 A 3 to the Hamiltonian in order to satisfy the equa- 
tions of motion for the new variables ip 2 and (p' 3 when tp' 2 = y? 3 = at time 0. 
The linear canonical transformation (A',ip') = (MA,M~ 1 cp) with 



M 



1 

1 



maps Hamiltonian (26) into Hamiltonian (24). 



3. 1 Torus with frequency vector oj 



We are looking at the break-up of the invariant torus with frequency vector u>o 
which is located at Q ■ A = for Hamiltonian (24) with e — 0. We notice that 
this torus is located at p = for Hamiltonian (23) and its frequency is equal to 
— 1. The three main resonances (of order e) are located at p — 1, p — —a, and 
p = — a 2 , and thus the torus is located in between the resonances p = —a and 
p — 1. We first discuss two rough estimates of the critical coupling obtained by 
some drastic simplifications : (a) Applying Chirikov's criterion [20] gives the 
following approximate value for the critical coupling e c ~ ((T+l) 2 /16 ~ 0.3377. 
(b) If we neglect the effect of the resonance located at p — —a 2 , e.g., by setting 
fj, — in the perturbation of Eq. (23), we can apply the renormalization 
procedure described in Ref. [21] for Hamiltonians with two degrees of freedom. 
Since Hamiltonian (26) does not depend on tp' 3 in this case, A 3 (t) is constant 
and the problem reduces to the study of the break-up of the invariant torus 
with frequency vector u>^ = (—1, er + 1) for the following Hamiltonian system 
with two degrees of freedom : 

H = ±A[ 2 + ■ (A[, A> 2 ) + £ (cos^ + costo + V ' 2 )). 
This method gives e c w 0.112. 



We perform the complete renormalization procedure described in the pre- 
ceding sections for Hamiltonian system (24) for a given value of /i. We fix 
the cut-off parameters L and J of the renormalization transformation, and 
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we take successively larger couplings e in order to determine whether the 
renormalization transformation 1Z converges to an integrable Hamiltonian, or 
whether it diverges. By a bisection procedure, we determine the critical cou- 
pling e c (/i] L, J) such that as n tends to +00 : 

7Z n H £ -> H for |e| < e c , 
TZ n H £ — > 00 for |e| > e c . 

Table 1 gives the values of e c for \x = 1 as a function of L and J com- 
puted by the two renormalization transformations. We notice that the values 
e c (/i = 1; L, J) converge to e c ~ 0.0886 as L and J grow. 
Figure 1 shows the values of e c (fi) for \x G [0, 10] determined by the renor- 
malization transformation for quadratic Hamiltonians with cut-off parameter 
L = 3. 

3.2 Frequency Map Analysis 

Frequency Map Analysis [5-8] allows to study the destruction of KAM tori by 
looking at the regularity of the frequency map defined from the action like 
variables to the numerically determined frequencies. In the present case (23), 
the frequency map is very simple as the system is equivalent to a one degree 
of freedom system with a quasiperiodic perturbation [8]. The angle variable 
x can then be fixed to x — 0, and we are left to a one dimensional map 
F : R — ► R, p — > v, where v is determined numerically from the output 
of the numerical integration of (p(t),x(t)) over a time interval of length T, 
starting with initial conditions (p(0),x(0)) = (po,0) [8]. On the set of KAM 
tori F is regular, or more precisely, it can be extended to a smooth function. 
Thus, when F appears to be non regular, this is an indication that all tori are 
destroyed in the corresponding interval. This allows to obtain a global vision 
of the dynamics of the system, as illustrated by figure 2 for Hamiltonian (23) 
with fi = 1, where v = F(p ) is plotted versus po for e = 0.09 and e = 0.12, 
and for p G [—2.5, 2.5] for a moderate precision (T = 1000). It seems clear on 
these figures that in the region between the two resonances v = — 1 — a 2 and 
v = 0, there are no tori left for e = 0.12, while many of them remain in the 
region v G [-1, 0] for e = 0.09. 

In order to have a more detailed view for the destruction of the tori with 
frequency v = —1, we have extended the time interval to T = 500000, and 
reduced very much the stepsize in p (figure 3). With these settings and for 
/j, = 1, it appears clearly that the torus with frequency v — — 1 is destroyed for 
e = 0.0895 (b), while the behavior of the frequency map appears to be very 
regular for e = 0.089 (a). It should be noted that the frequency map analysis 
provides a criterion for the destruction of tori. The fact that the frequency 
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curve appears to be irregular provides an evidence that the tori are destroyed, 
but when the curve is regular, a higher accuracy (which means a longer time 
span T) could reveal the destruction of additional tori. 

For e = 0.0895 and for /j, — 1, all tori in the vicinity of v — — 1 are destroyed, 
which is in agreement with the value e c ~ 0.0886 found with the renormaliza- 
tion technique, as this is the largest value of e for which the renormalization 
converges. 

We have performed computations for other values of the parameter /i. Fig- 
ure 1 shows the agreement between the couplings obtained by Frequency Map 
Analysis and the ones obtained by the renormalization transformation. A bet- 
ter accuracy can be obtained by taking a larger cut-off parameter L for the 
renormalization computations. 



3.3 Last K AM torus 



It appears clearly in figure 2 that the torus with frequency v — — 1 is not 
the last torus to be destroyed in the interval of frequencies [—1 — <r 2 ,0] for 
/i = 1, and that many tori still survive for e > 0.0895. We thus have searched 
for the value of the parameter e for which the last torus disappears. This is 
done by increasing the value of the parameter e until all tori disappear. In this 
procedure, at each stage we identify the resonant islands and chaotic regions, 
and then decrease the stepsize in p . This allows one to obtain very refined 
details as one gets close to the critical value of e. 

The frequencies of the last invariant torus can be investigated by frequency 
map analysis. This method gives a critical threshold at about e c ~ 0.11 (see 
figure 4). For Hamiltonian (23) with fx — l, the last invariant torus is located 
between the resonances p — 1 and p = —a at p ss 0.354. The corresponding 
frequencies for the last three-dimensional torus for the family of Hamiltoni- 
ans (24) are po + a 2 , po + a and 1 — p . 

A nearby invariant torus has the frequencies (a + 2)(l — p ), (a' 2 + 2)(1 — p ) 
and 1 — Po- The critical coupling for the break-up of this invariant torus can 
be computed by the renormalization transformation defined in Sec. 2, mainly 
by defining first a unimodular transformation that maps the frequencies of the 
torus into the frequencies a 2 , a and 1 (the perturbation is then different from 
Eq. (25)). The renormalization gives e c m 0.11. 

If we neglect the resonance located at p — —a 2 , renormalization methods and 
frequency map analysis show that there are two last KAM tori, one located 
at p ~ 0.112 and the other one at p ~ —0.436, with a critical threshold of 
e c ~ 0.149. Thus the effect of the third resonance at p — —a 2 is to destabilize 
the motion in the region between p — 1 and p = —a closer to p — —a (since 
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the third resonance is closer to p — — er). Consequently the last KAM torus 
is located nearer p — 1 than in the system without the third resonance, and 
the value of the threshold of global stochasticity (break-up of the last KAM 
surface) is smaller. 
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Table 1 

Critical coupling e c for Hamiltonian (23) with n = 1 as a function of the cut-off 
parameters L and J, computed with the transformation RGl defined in Sec. 2.1, and 
with the transformation RG2 defined in Sec. 2.2. Frequency Map Analysis indicates 
that the spiral mean torus is broken for couplings larger than 0.0895. 
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Fig. 1. Values of the critical couplings £ c (/-0 for Hamiltonian (23) with fi G [0,10]. 
The dashed curve is obtained by the renormalization method, and the continuous 
curve is obtained by Frequency Map Analysis. 
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Fig. 2. Frequency map for Hamiltonian (23) with jjl = 1 and with e = 0.09 (a) and 
e = 0.12 (b). The dotted lines correspond to v = 0,-1,-1 - a, -I - a 2 . T = 1000. 
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Fig. 3. Frequency map for Hamiltonian (23) with [i = 1 and with parameter 
e = 0.089 (a) and e = 0.0895 (b). v + 1 is plotted versus (po — poo) x 10 10 , where 
Poo = -3.3307 x 10~ 4 . The dotted lines corresponds to v = -1. T = 500000. 
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Fig. 4. Determination of the last KAM tori for the Hamiltonian (23) with fj, = 1. 
For each selected value of the parameter e, the set of frequencies v for which tori 
are not destroyed is plotted versus e. The critical value of e for which no tori will 
survive is thus e c « 0.11. 
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